Under-the-barrier dynamics in laser-induced relativistic tunneling 
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The tunneling dynamics in relativistic strong-field ionization is investigated with the aim to develop an intuitive 
picture for the relativistic tunneling regime. We demonstrate that the tunneling picture applies also in the 
relativistic regime by introducing position dependent energy levels. The quantum dynamics in the classically 
forbidden region features two time scales, the typical time that characterizes the probability density's decay of 
the ionizing electron under the barrier (Keldysh time) and the time interval which the electron spends inside 
the barrier (Eisenbud-Wigner-Smith tunneling time). In the relativistic regime, an electron momentum shift as 
well as a spatial shift along the laser propagation direction arise during the under-the-barrier motion which are 
caused by the laser magnetic field induced Lorentz force. The momentum shift is proportional to the Keldysh 
time, while the wave-packet's spatial drift is proportional to the Eisenbud-Wigner-Smith time. The signature of 
the momentum shift is shown to be present in the ionization spectrum at the detector and, therefore, observable 
experimentally. In contrast, the signature of the Eisenbud-Wigner-Smith time delay disappears at far distances 
for pure quasistatic tunneling dynamics. 
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Introduction With present-day strong lasers providing in- 
tensities of the order of 10 22 W/cm 2 (T| and employing highly 
charged ions, the relativistic regime of strong field atomic pro- 
cesses, in particular, strong field ionization of highly charged 
ions 0, is accessible Q. Relativistic generalizations of the 
strong field approximation (SFA) (3H3 and the quasiclassical 
theory of ionization [9) provide simple ways to calculate 
ionization probabilities in the relativistic regime ifTOlfTTTl . 

The ionization occurs in the so-called tunneling regime 
when the laser's electric field amplitude E Q is smaller than those 
required for over-the-barrier ionization lfT2l and the Keldysh 
parameter y = oj *^2T p jE§ lf]~3l is well below unity with the 
ionization potential I p and the laser's angular frequency oj. 
The Keldysh parameter defines the time scale tk = y/oj during 
which individual momentum components of the ionized wave 
packet are formed |5|. When this time scale is shorter than 
the laser period, i. e., y <?c 1, the laser field can be treated as 
quasistatic. 

In the nonrelativistic as well as in the relativistic case, the 
electron travels during its journey from the bound state into the 
continuum through a region which is forbidden according to 
classical mechanics. For nonrelativistic dynamics, when the 
typical electron velocities are negligible compared to the speed 
of light, it is legitimate to treat the laser field as a homogeneous 
time-dependent electric field E(t) and to neglect its magnetic 
component. In this case a well-known intuitive picture for 
the tunneling regime arises in which the electron tunnels out 
through the effective potential V e ff (T, = V(r) + r • E(t) formed 
by the atomic potential V(r) and the laser field. The Keldysh 
time tk may be interpreted as the time that a classical free 
electron would need to cover the length / ~ I p /E of the tun- 
neling barrier at the characteristic velocity of a bound electron 
v a = ^2T p . The contribution of the laser's magnetic field com- 
ponent becomes non-negligible when we enter the relativistic 
regime. A transverse laser field with perpendicular electric and 
magnetic components cannot be described solely by a scalar 



potential. Consequently, a generalization of the tunneling pic- 
ture into the relativistic regime is not straightforward ifTH and 
there is no clear intuitive picture for the relativistic quasistatic 
ionization dynamics in the classically forbidden region. 

Later, we will identify the so-called tunneling time r E as a 
second time scale that is fundamental for the electron dynamics 
in the classically forbidden region. The tunneling time through 
a potential barrier (Eisenbud-Wigner-Smith time delay (T5J) 
has been a long-standing problem in quantum mechanics ITT51 — 
[T8l . Interest to this problem has been renewed recently by 
experimental attempts to measure the temporal delay at the 
electron emergence into the continuum from a bound state 
during the laser-induced (nonrelativistic) tunneling |fT9l . A 
vanishing tunneling time within the experimental error has 
been observed in these experiments. 

In order to develop an intuitive picture for the relativistic 
tunneling regime of strong field ionization, we carry out a 
detailed analysis of the electron dynamics through the classi- 
cally forbidden region in the relativistic regime and compare it 
with the corresponding nonrelativistic case. The characteristic 
time scales are examined. Our analysis is based on the SFA 
EH3 as well as on the exact solution of the relativistic wave 
equations and numerical simulations. The possibility for the 
experimental detection of the specific relativistic features of 
the classically forbidden dynamics is also elaborated. 

System A highly charged hydrogenlike ion with nuclear 
charge Z < 100 at the origin of our coordinate system is 
considered to interact with a strong linearly polarized laser 
field with electric field of amplitude Eq in x-, magnetic field 
of amplitude Bq in y- and propagation in z-direction (Eo = Bo). 
Position and momenta of electrons that are ionized by the 
laser pulse may be measured by a detector far away from the 
interaction zone. Throughout this Letter we describe the laser 
field via its vector potential in electric field gauge |20|[2T1, i. e., 
0(r, t) = -r • E(r]) and A(r, t) = -k[r • E(t])]/c with the laser 
phase r] = t - z/c, the electron's coordinate r = (x,y,z), the 
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unit vector in laser propagation direction k, and the speed of 
light c. The dynamics of the electron state |^) in the rest frame 
of the ionic core is then defined via the Dirac equation 

id t |¥> = [ca • \p + Air, t)] + V(r) - 0(r, t) + /3c 2 ] |Y> (1) 

with the momentum operator /) and the Dirac matrices a and 
J3. In the following, we use an extreme but feasible set of 
parameters, viz. I p /c 2 = 0.25 and Eo/E a = 1/30, with E a = Z 3 
denoting the characteristic atomic field strength, representing a 
hydrogenlike ion with nuclear charge of Z « 90 in the tunneling 
regime (T2) . 

Relativistic strong field approximation The SFA al- 
lows us to calculate the ionized part of the asymptotic electron 
wave packet at a detector in the remote future. As we are in- 
terested in the under- the-barrier dynamics of the electron, the 
momentum distribution at the tunnel exit time to is deduced 
analytically from the distribution at the detector, assuming free 
propagation in the laser field outside the barrier. The asymp- 
totic ionization amplitude in momentum space reads (22) 



(X 



(i/s Pf , Sf \r 'E(rj)\{fr ) dq , 
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where pf = (Pf, x ,Pf, y ,Pf, z ) and Sf denote the electron's final 
momentum and spin, respectively, \ifs p , s ) is the Volkov state 
(23) with momentum p and spin s, while \\J/q) is the dressed 
Dirac bound state spinor of a short range potential with binding 
energy -I p (22). The relativistic wave function \i//q) is used 
in the relativistic as well as in the nonrelativistic calculation 
allowing to trace back the features in the under-the-barrier dy- 
namics to the Hamiltonian's different corrections to the electric 
dipole approximation. Calculating the integral ^ via the sad- 
dle point method, the momentum distributions at the detector, 
see Fig.[T](a), and at the tunnel exit are derived; see Fig.[T](b). 
While parabolic wings of size pf fZ ~ E 2 /(qj 2 c) are a common 
feature of asymptotic momentum distributions for ionization 
in the relativistic regime (22), Fig. [T](a) reveals also a shift 
of the global maximum of the asymptotic momentum distri- 
bution from (pf, x >Pf,z) = (0,0) in the nonrelativistic regime 
to (0, I p /(3c)) in the relativistic regime (24) . The shift of the 
global maximum arises during the under-the-barrier motion. 
This may be realized by calculating the distribution of the ki- 
netic momentum at the exit nonrelativistically and comparing 
it with the fully relativistic result; see Fig.[T](b). To identify 
the relevant relativistic effects, we expand the Dirac equation 
([]]) retaining the magnetic dipole correction and the leading 
relativistic correction to the electron's kinetic energy (termed 
mass -shift correction) 
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+ V(r) + r • E(t) \ |Y> . 



(3) 



This expansion is justified because relativistic effects for the 
under-the-barrier motion are determined by the parameter I p /c 2 



and even for the extreme parameter set used the tunneling 
process is only weakly relativistic. This can be seen from 
Fig.[T](b), where the calculation based on Eq. ^ recovers the 
exact result from Eq. ([T]). Furthermore, the figure shows that 
the magnetic dipole term (Lorentz-force) is responsible for the 
momentum shift and a decrease of the tunneling probability 
whereas the mass-shift correction counteracts the latter effect. 

A quasiclassical picture for relativistic tunneling 
For a qualitative interpretation of the tunneling procedure 
(y « 1) and the momentum shift, we omit the mass-shift 
correction as justified above and employ the quasistatic ap- 
proximation in the quasiclassical picture. Then the Hamilton 
function corresponding to Eq. ^ yields at maximal laser in- 
tensity 



Px Py 
2 2 



(p z + xE /c) 2 



+ V(r) - xE , 



(4) 



with the canonical momenta p x , p y , and p z , r = ^x 2 + y 2 + z 2 . 
In ^ the electric field gauge [20 , 21 ] ensures that the Hamilton 
function equals the total energy -I p . Similarly as shown in 
the nonrelativistic case 1 12], we estimated numerically that 
the relativistic deviations still result in tunneling dynamics in 
close vicinity along the x-axis (the polarization direction). The 
atomic potential's central force in the direction of the j-axis 
and the z-axis (the laser's propagation directions) is then also 
negligible in the tunneling region and the canonical momenta 
p y and p z become approximately conserved. Consequently, the 
electron dynamics is separable and the equation H = -I p can 
be written as 



Py ( Pz +XE /C) 2 

2 2 



= y + [V(x)-xE ] 



(5) 



Equation §5§ is a generalization of an often applied nonrela- 
tivistic one-dimensional model (12). It may be interpreted as 
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FIG. 1 : (a) Density plot of the asymptotic spin-averaged electron 
momentum distribution at the detector. The distribution's maximum 
(white cross) is shifted from p fz = (nonrelativistic limit, white line) 
to pf jZ - I p /(3c). (b) The distribution of the electron's kinetic mo- 
mentum in the laser propagation direction at the tunnel exit at time t , 
with maximal laser electric field, in the electric dipole approximation 
(green, dotted line), including the magnetic dipole correction (blue, 
dashed line) and the magnetic dipole plus the mass- shift correction 
(red, dot-dashed line) and fully relativistically via the Dirac equation 
in the SFA (black, solid line). The applied parameters are I p /c 2 = 0.25 
and E /E a = 1/30. 
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FIG. 2: (color online) Scheme of laser induced tunneling in a 
Coulomb potential with I p /c 2 = 0.25. The potential barrier V(x)-xE 
(black line) is displayed with E /E a = 1/30 and the relativistic and 
nonrelativistic modified energy levels. The dotted green line indicates 
the description in electric dipole approximation, whereas magnetic 
dipole effects are taken into account for the solid red line. The mo- 
menta p y and p z are chosen such that the WKB tunneling probability 
is maximal. 
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FIG. 3: (a) The scaled time delay vs. the coordinate along the laser 
polarization direction, via a quantum mechanical (black, solid line) 
and a quasiclassical description (gray, dashed line), (b) The electron 
wave packet trajectory in the x-z-plane, via a quantum mechanical 
(black, solid line) and a quasiclassical description (gray, dashed line). 
In both subfigures x/x e < 1 corresponds to under- the-barrier motion 
with x e = Ip/Eo denoting the tunnel exit coordinate. The applied 
parameters are I p /c 2 = 0.25, E /E a = 1/30, and the most probable 
initial kinetic momenta p z ,km(xd = -2I p /(3c) and p y = 0. 



follows. During the quasi one-dimensional motion along the 
laser polarization direction, the electron propagates through 
the barrier V e ff(x) = V(x) - xEq with a position dependent 
total energy s(x) = -I p - {p 2 /2 + [p z + xEq/c] 2 /2], which 
is the difference between the binding energy and the kinetic 
energy of the motion in the propagation direction of the laser. 
Equation ^ yields the momentum in electric field direction 
p x = i V2[Veff(x) - s(x)]. Note, that the kinetic momentum in 
the propagation direction of the laser p z ^(x) = p z + xEq/c in 
the laser's propagation direction is shifted due to the Lorentz 
force during the under-the-barrier motion causing a modifica- 
tion of s(x) and, consequently, a modification of the momentum 
p x as compared to the nonrelativistic case where the kinetic 
momentum p Zf \dn(x) = Pz is conserved. The WKB tunneling 
probability is proportional to exp(-r), where T is given by 
the integral Y = -2i f e p x dx over the classical forbidden re- 
gion Xi< x < x e El [Hi. For short range potentials the WKB 
tunneling probability is maximal when tunneling starts with 
a negative kinetic momentum p z ^ n (xi) - -2I p /(3c) reaching 
the tunnel exit with positive p z ^m(x e ) = I p /(3c). At this opti- 
mal momentum the energy s(x) is strictly below -I p , which is 
the corresponding expression in the case of the electric dipole 
approximation, see Fig. [2] This indicates a decrease of the 
maximal tunnel probability with respect to the nonrelativistic 
case, cf. with Fig.[T](b). 

In the WKB approximation the electron wave function 
evolves along quasiclassical trajectories JHJ |TT||. Under the 
barrier, they are described by classical Hamiltonian mechanics 
developing in imaginary time. The imaginary time determines 
the exponential damping of the ionized electron wave function 
during tunneling. In particular, a time interval -\tk is required 
to pass the barrier; thus, the wave function is reduced by a 
factor e"^ T/: while tunneling. Under the barrier the Lorentz 
force of order v a Bo/c acting for the Keldysh time tk causes 
a total momentum shift of A/? z?k i n = p z ^ n (x e ) - p z ^ n (Xi) ~ 
TKV a Bo/c ~ (x e - Xi)Eo/c ~ I p /c. In order to compensate for 
the acceleration in positive z-direction by the Lorentz force, the 



z-component of the initial kinetic momentum at entering the 
barrier has to be negative similar to l25l . The kinetic momen- 
tum in the propagation direction of the laser is positive at the 
barrier exit, see the dashed blue lines in Fig.[T](b), in contrast 
to the nonrelativistic case where it vanishes. Considering the 
mass-shift correction in addition to the magnetic dipole terms 
increases the tunneling probability; see Fig.[T](b). The general 
tunneling picture, however, is not affected. 

Tunneling time We point out that although the Lorentz 
force induces a momentum shift along the laser propagation 
direction during the under-the-barrier motion, nevertheless, in 
the quasistatic WKB approximation the tunneling process hap- 
pens instantaneously and thus there is no real coordinate drift in 
laser propagation direction as a result of the under-the-barrier 
motion. However, in reality the under-the-barrier motion re- 
quires a small nonzero time span which can be derived by 
going beyond the quasistatic WKB theory. 

For this purpose we use the exact wave function of the elec- 
tron in the quasistatic potential formed by the laser field and 
a zero-range atomic potential. The kinetic energy term of the 
Schrodinger equation with the magnetic dipole corrections is 
expanded in the coordinate x around the tunnel exit coordi- 
nate x e . Then, the electron wave function in the continuum is 
derived as 

i/f + (x) = TAi {e- /3 [2 1/3 (xE - I p ) /Ef + cf>(x)]} , (6) 

where <j>(x) = (p z c + I p )(3p z c + 2xE + I p )/3c 2 (2E ) 2/3 , 
and Ai denotes the Airy function. The real trans- 
mission coefficient T is obtained by matching if/+(x) 
with the bound state wave function fTTTl i//q(x) = 
VSnWfrfl + (Py + P 2 z )/(2v a )] exp {-(v a + (p 2 y + p 2 z )/(2v a )]x}. 

The time it takes for the maximum of the wave packet to 
reach position x [Eisenbud-Wigner- Smith time delay r E (x)] 
and the drift of the electron wave packet in propagation di- 
rection at position x during the under-the-barrier motion are 
determined by the phase of the wave function arg[^ + (x)] fTTlL 
viz. r E (x) = -darg[ifr+(x)]/dl p andz(x) = -darg[ifr+(x)]/dp z . 
Figure [3] (a) shows that the electron spends some time un- 
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der the barrier when it is close to the exit. The tunneling 
time, i. e., the time delay at the tunneling exit is given by 
r E (x e ) * 3 5/6 ^(E a /E ) 2/3 /[2T(l/6)I p ] and the drift distance 
at the tunneling exit is z(x e ) * ^/ p /[3 1/6 r(l/6)^Q /3 c]. The 
tunneling time of order 1 /I p may also be estimated via Heisen- 
berg's uncertainty relation when assuming an energy uncer- 
tainty of order of I p . The drift size along the laser propagation 
direction at the tunneling exit is proportional to the tunneling 
time with a proportionality factor that is equal to the kinetic mo- 
mentum in the laser's propagation direction at the tunnel exit, 
z(x e ) = p z ,\dn( x e)TE(Xe)- This indicates that the drift is shaped 
during the Eisenbud-Wigner- Smith tunneling time. Note that 
for the applied parameters the drift size at the tunneling exit 
z(x e ) is of order of 1 jc. The drift along the laser propaga- 
tion direction during the under-the-barrier motion is due to the 
Lorentz force and disappears in the nonrelativistic limit. The 
time delay r E (x), however, behaves similarly to that shown in 
Fig. [3] (a) also in the nonrelativistic regime. 

The time delay far from the tunnel exit, that means x - x e > 
£"~ 1/3 , follows approximately r E (x)\ x ^oo ~ ^2(x - x c )/Eq. 
The latter corresponds to the quasiclassical motion time of 
the electron which tunnels instantaneously and continues at x e 
with vanishing initial momentum in the laser's electric field 
direction p x = 0, reaching x at a moment t (x - x e = Erf 2 / 2). 
Therefore, we can conclude that the signature of the nonzero 
tunneling time exists only in a small area near the tunnel exit 
at distances Sx ~ E~ l/3 and disappears at far distances near the 
detector, see Fig.[3](b). This consequence is relevant also in 
the nonrelativistic case and explains the observation of the van- 
ishing tunneling time in the experiments lfT9l . Mathematically, 
the problem of the vanishing Wigner time at large distances 
arises when the laser-induced quasistatic barrier is linear in the 
x coordinate around the exit of the tunneling electron in the in- 
terval [x e - Sx, x e + Sx], which is the case when the saddle point 
of the effective potential barrier formed by the laser and atomic 
potentials is far from the tunneling exit. Instead, in the regime 
close to over-the-barrier ionization at Eo/E a > 1/25 the effec- 
tive potential has a parabolic character and analog calculations 
as carried out here indicate a measurable Eisenbud-Wigner- 
Smith tunneling time l26l . In interpreting experimental results 
on this tunneling time, the distortion of the electron trajectories 
by the Coulomb field of the ionic core during the sub-barrier as 
well as out-of-the-barrier motion should be properly accounted 
for, see Ref. 1271 . 

Numerical simulations To corroborate our analytical 
arguments, we have carried out an ab initio simulation of the 
tunneling process in a highly charged ion in a laser field of 
relativistic intensities by solving the time-dependent Dirac 
Eq. ([T]) numerically (28). Figure [4] shows the wave packet 
of the active electron at the moment when the laser electric 
field is at its maximum. Here one part of the wave packet is 
still bound or under the barrier while the other part is already 
in the continuum. The part under the barrier is not exactly 
symmetric to the laser polarization direction. There is a small 
drift motion into the laser propagation direction which is of the 
order z ~ 1 /c in accordance with our estimation above. 
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FIG. 4: (color online) Numerical simulation of the electron wave 
function in a soft-core potential via the Dirac equation. The density 
plot shows the electron density at the moment when maximal laser 
field strength is attained. The solid black line indicates the maximum 
of the density in the laser propagation direction while the dashed 
line corresponds to the most probable trajectory resulting from the 
quasiclassical description. Solid green lines correspond to the border 
of the classical forbidden region. White arrows and the cross indicate 
the directions of the laser's electromagnetic fields and its propagation 
direction. The inset shows a scale-up of the region close to the tunnel 
exit. The applied parameters are I p /c 2 = 0.25, and E /E a = 1/30, 
and y = 0.1. 



Conclusion Laser-driven relativistic tunneling dynamics 
was shown to be understood by a simple picture incorporat- 
ing a scalar tunneling barrier and position-dependent energy 
levels. A spatial shift and a momentum shift along the laser 
propagation direction at the tunnel exit were identified as signa- 
tures of the relativistic dynamics through the tunneling barrier. 
The spatial drift is shown to be proportional to the Eisenbud- 
Wigner- Smith time delay, while the momentum shift to the 
Keldysh time. While the spatial shift features intrinsic limita- 
tions for its measurement, the momentum shift is of the order 
of 10 a. u. for the applied parameters and within the reach of 
up-to-date experimental techniques fT9l . 

We are grateful to C. Muller for valuable discussions at the 
onset of the project. 
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